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Abstract 

We investigate dynamical systems with time-dependent mass and fre- 
quency, with particular attention on models attaining the minimum value 
of uncertainty formula. A criterium of minimum uncertainty is presented 
and illustrated by means of explicit and exactly solved examples. The role 
of the Bogolubov coefficients, in general and in the context of minimum 
uncertainty case, is discussed. 

PACS: 03.65.-w; 03.65.Ca; 42.50.Dv 

1 Introduction 

We analyze some peculiar situations emerging in the study of oscillators charac- 
terized by constant and time-dependent mass and frequencies. Our investigation 
is focused mainly on the treatment of a nonlinear auxiliary equation which can 
be associated with the systems under consideration ( 11161 1. A special attention 
is paid on the problem of the formulation of the uncertainty relation, where in 
our framework the auxiliary equation has a basic role. 

A classical generalized or time-dependent oscillator (TDO) is governed by 
the Hamiltonian 

H{t) = ^ + \mu:W (1.1) 

where q,p are conjugate variables, and m — ni{t), uj = Lo{t) are given functions 
of time. The Hamiltonian Hl.l|l gives rise to the equation of motion 

q + Mq + <jj'^q = Q , (1.2) 
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where M = M{t) — rh / m and dot means time derivative. Through the trans- 
formation 

-i f* M(t')dt' o\ 

e 2 -"o ^ ' y ^ (1.3) 

Equation p.2f) can be cast into the form 

y + n^{t)y = 0, (1.4) 



where 



n^{t) = - (auj^ - 2M - M^) . (1.5) 



4 



The foUowing equation 



K 



a + n\t)a:^ — , (1.6) 

where X is a constant, can be related to Eq. (|1.4|) . in the sense that if yi and 
2/2 are two independent solutions of Eq. (|1.4|) , then the general solution of the 
auxiliary equation (|1.6|) can be written as 0] 

a^iAyl + Byl + 2Cym)K (1-7) 
A, B, C being constants such that 

AB^C^ = ^ (1.8) 

and Wo = 2/12)2 — 2/i2/2 = const is the Wronskian. From the theory of the 
auxiliary equation (|l.t)|) a phase is involved given by a real function 6{t) defined 

by 



m = / ■ (1-9) 



dt' 
It. ^) 

(See jll3l7| : for some applications: 07 El)- 

In this Letter we discuss some aspects of the quantum theory of generalized 
oscillators with time-dependent mass and frequency. Precisely, we revise the 
formulation of the uncertainty relation in terms of solutions of the auxiliary 
equation H1.6I) . with a particular care to the analysis of the Bogolubov coeffi- 
cients and the possibility to attain the minimum value of the uncertainty product 
(AqQ) (Aq,P), where AqQ and A^P denote the variances of the position and 
momentum operators between extended coherent states. A simple critcrium is 
found which allows one to build up time-dependent models reaching the mini- 
mum uncertainty relation. Some examples of these systems are presented. One 
of these, which leads to an equation of the Bessel-type, is associated with an 
auxiliary equation (see H4.10|l 'l exactly solved in terms of a (convergent) power 
series expansion. 

The outline is as follows. In Section 2 some examples of (real) exact solu- 
tions of the auxiliary equation H1.6|l are shown. In Section 3, the uncertainty 
relation in terms of the Bogolubov coefficients of the generalized oscillator with 
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time-dependent mass and frequency is considered. A criterium of minimum un- 
certainty relation is also presented. Section 4 concerns with examples of exactly 
solvable TDO's minimizing the uncertainty formula. In Section 5 a concluding 
discussion is reported. 

2 Examples of exact solutions of Eq. ( 11.61 ) 

Two interesting cases where Eq. can be exactly solved via (|1.7|l . are 

represented by a) the usual (time-independent) harmonic oscillator and b) the 
Kanai-Caldirola oscillator. In what follows, we shall deal with both cases. 

Case a) The harmonic oscillator 

Let us assume that uj{t) and rn(t) are constants, namely uj{t) = ujq and 
m{t) = mo. Then, the Hamiltonian Hl.l|) takes the form 

Hit)^£- + lmo.U^ , (2.1) 

while the auxiliary equation becomes 

a + Lula^ . (2.2) 

The general solution of this equation can be obtained via (|1.7|l by choosing, say, 
Qi = 90 coso^o^, <Z2 = 9o sincjoi • If we put K = j for convenience, and hereafter 
we limit ourselves to the coefficients A,B,C leading to real solution (7{t), then 
equation (|1.8|l is satisfied by 

A = B = —, C = (2.3) 

for k — luq {k is defined below) , and 



B = ^-272 + ^o,,2„2 ° cos(2ci ± 2^oto) , (2.5) 



C = T ^n\ 2'^ ^sin(2ci±2a;o^o) , (2.6) 

where ci is an arbitrary constant, for k > loq, where A; is a constant of integration 
appearing in the equation 

' -^'^ ' ^ = k (2.7) 
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arising from Inserting the coefficients (|2.3|) and (|2.4|l - (|2.6|l into (|1.7|l 

provides 

cj - ^= , (2.8) 



and 



2^0 



2^0 



-2- cos(2ci ± 2u;oi) 



(2.9) 



respectively. The phases ()1.9|l corresponding to the solutions H2.8|l and l|2.9|l are 
given by 

9{t)^2ujo (t-to) 

and 



e{t) = ±2tan"^ 



UJQ 



tan (ci ± ujQt') 



(2.10) 
(2.11) 



respectively. 

Case b) The Kanai-Caldirola oscillator 

Now let us suppose that the frequency U'^{t) defined in Ijl.SII is a constant, 
i.e. fl^{t) — fig, and m(t) = uiq exp(7t) , where 7 = Mq is the damping. So the 
equation of motion H1.2|l corresponds to a particular form of the Kanai-Caldirola 
Hamiltonian |1UI11| . i.e. 

„2 1 



H 



exp (-7 1) + ^mo exp (7 1) fll 



The auxiliary equation l|1.6|l can be written as 

1 



4a3 



Since from (|1.5|) we get 



(2.12) 



(2.13) 



(2.14) 



we have two possible situations: fig > and ii) Vlf^ < Q . In case i) two 
independent solutions of the equation of motion (|1.2|) are qi — q^ cos fipt , 
92 = 9o sin rigi ■ Consequently, case i) is similar to the case A/ = considered 
previously. Therefore, for a we find two possible solutions, which coincide with 
H2.8|l and H2.9|l where is replaced by fio- Conversely, case ii) is concerned 
with the solution 



ci 



which implies the phase 



e{t) = ±2tan"i 



1 



4irjn 



cosh(2c2±2|f7o| t) 



(2.15) 



4 |f7o|^ c2 + 1 - 2 ir^ol ci I tanh (c2 ± jr^ol 



(2.16) 
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3 The uncertainty relation in terms of the Bo- 
golubov coefficients for the TDO 

In 1^, in the context of the determination of Noether invariant operators and 
dynamical group of the TDO described by the Haniiltonian operator 



(3.1) 



the foUowing time-dependent lowering and raising operators a{t) and a^(i) have 
been introduced: 



TON 1 


r Q 




[2a 




r Q 


li) 


I2a 



m \ 2 



a 



-P- a 



M 



a\Q 



a{t) = 
a\t) = 

\ It / I Z,C I III 

where m = m{t), M ^ ^, [Q, P] = ih, [a{t),a^t)] = 1, and 

1 



4cr3 



where $7^ is given by H1.5I) . Now, from (|3.4|l we obtain 



where F is defined by 

F = J niia^dt' = ( 

Form 1(221 and (E^J) we find 

Q = \ —<j{a + a^) , 
V m 



4:Louj - M - M M] dt' 



P = Vhm {(a + ^ 



where 



and 



i 








.2^ + 





fM V' 








— (T — a 




k + F + 




V 2 J 







Now, let us deal with the variances 

(A„Qf = (g2)-(g)2 , {A^pf = {p')-{pf 



(3.2) 
(3.3) 

(3.4) 
(3.5) 

(3.6) 

(3.7) 
(3.8) 

(3.9) 

(3.10) 
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where the expectation value is given by ( 
generalized coherent state. Thus, we have 



I a) denoting the 
(3.11) 



which entail the uncertainty formula 



(A„Q) (A„P) = - 



ha 



M 

1 + 4(T^ ( cr - —a 



k + F- 



M 



.f2 



n^]a^ - Man 



(3.12) 



which holds for a generalized oscillator with time-dependent mass and frequency. 
For m = mo {i.e. M = 0), Eq. 1)3.12(1 reproduces the relation (see 



(A„Q) (A„P) = -v/l + 4a2a2 > - , 



(3.13) 



where the minimum is attained only for tr = 0. 

In order to write the Bogolubov transformation, we introduce the operators 



Q 



. hmoLUQ ( I 
P = — 7, — [ao- % 



(3.14) 



where oq = ao(0) is the Schrodinger picture fixed photon annihilator operator. 
Substitution from 1)3. 14() into 1)3. 2|l and ()3.3)l provides 



a{t) ^ ^i{t)aQ + v{t)al , 
where the Bogolubov coefficients /i(t) and v{t) take the form 



^i{t) 

with 



moLUo ■ 

7] H CT 



2moWo 



rriQUJo 
•q (7 



1 -f^i 



For any real solution of Eq. 1)3.4(1 we get 



k + F- 



k + F- 



m? 2 



4^ I cr^ - Merer 



'0^0 



M2 



M 



cr"' - Maa 



(3.15) 

(3.16) 
(3.17) 



1 

^2 
(3.18) 
1 

~ 2 
(3.19) 
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so that the complex functions ^(i) and ^{t) possess the property 

Ht)f-W{t)\' = l . (3.20) 

Furthermore, the quantities 



,it) + .(t) . , f4t) .it) = J'-^a (3.21) 

V moUJo V m 

can be involved in the uncertainty formula (|3.12|) , i.e. 



(A„Q) (A„P) ^ ^ 



1 + 4(7^ f cr - ^(T 



= 2 (4^'h|')'=-|M(t) + K^)IKi)-Kt)l> 2 ■ (3-22) 

The uncertainty formula H3.22|l is closely related to the concept of coherent states 
for the generalized oscillators. We remind the reader that such coherent states 
in the context of Lewis-Riesenfeld theory were constructed by Hartley and Ray 
in 1982 These states share all the features of the coherent states of the 
conventional (time-independent) oscillator except that the uncertainty formula, 
i.e. the product of position and momentum is not minimum. A few years later, 
Pedrosa showed that the coherent states devised by Hartley and Ray for 
the TDO are actually equivalent to the well-known squeezed states (see, for 
instance, ^31 ^^^d llSj). 

A criterium of minimum uncertainty for a time-dependent oscillator with 
variable mass and frequency is expounded in next Subsection. 

3.1 A criterium of minimum uncertainty relation for the 
TDO 

The product l|3.12|l reaches its minimum value whenever the condition 

^a = a (3.23) 

is fulfilled. This implies the constraint 

a = , (3.24) 

where c is a constant of integration. By introducing (|3.24() in the auxiliary 
equation (|1.6(l yields 

2mm-m'^ + 4:n'^m'^ = \ . (3.25) 
This equation can be elaborated as follows. Substitution from 

mm^m^- h (3.26) 

dt m 
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into 13.25II . wc find 

m{t)u;{t) = ^ (3.27) 

by virtue of 1)1. 5|l . Equation (|3.27|) represents a criterium of minimum uncer- 
tainty for the product of variances (A^Q) (A^P) for the generahzed osciUator 
1)1. 2|l . That is, an oscillator with generally time-dependent mass and frequency 
minimizing the uncertainty formula is described by the equation 

q--q + u^q^O . (3.28) 

UJ 

In the light of this result, it is clear once again that for the usual harmonic 
oscillator (corresponding to the particular case in which oj{t) = = const 
and m = Too = (2woC^)^^), the unique exact solution of the auxiliary equation 
minimizing the uncertainty product is cr = . The other solution H2.9|l can 

minimize the uncertainty formula only approximately, say for ljq < k << 3oJq. 
Notice that in the minimum uncertainty case the Eq. H1.9|l takes the form 

9{t) = 2 f Lu{t')dt' . (3.29) 

J to 

4 Examples of exactly solvable TDO's minimiz- 
ing the uncertainty formula 

Below we shall display a few interesting examples of TDO's which correspond 
to the minimum of the uncertainty formula. 

Case I) A TDO with exponentially decreasing frequency 

Let us consider a TDO with frequency 

a; = Wo exp (-7oi) , (4.1) 

and mass 

^ exp(7ot) , (4.2) 



2c-^a;o 

so that the criterium (|3.27|) is verified (70 = const > 0). The equation of motion 
H1.2|l reads 

g + 7o9 + t^o exp (-270 t) q = . (4.3) 

By means of the transformation (|1.3|) . i.e.q — exp ("y^) Equation 14. 3|) can 
be cast into the equation 

y + n^t)y^O , n^t)^u;lexpi-2jot)-^ , (4.4) 
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as one finds from H1.5|l . The general solution of Equation H4.4|l can be given in 
terms of elementary functions: 



y = V7oe 



,70 t 



ci COS 1 — e-"^"* 
70 



C2sml ^e-^"* 
.70 



where Ci and C2 are arbitrary constants. 



Finally, due to Eq. (1223, the phase is = -2^e' 



70* 



Case II) A TDO with mass m = mot^ and frequency lo 



t 

to 

1 



In this case the equation of motion H1.2|l is 



1 



q + -q 



(2moc2t2)2 
which may be transformed into 

1 



9 = 



y 



■y = 



via q — J- Equation (|4.5f) allows therefore the general solution 



ci cos ■ 



C2 sm ■ 



(2moc2)t (2moc2)t_ 
Ci, C2 being arbitrary constants. From Eq. H3.29|l we also get 9 - 
Case III) A TDO described by an equation of the Bessel-type 



(4.5) 



(4.6) 



(4.7) 



-1 



Let us suppose that 



m = mo a{t) 



UJ 



^0 

a{t) 



Then, muj = motoo = (2c ) . If we demand that 



(4.8) 



(4.9) 



where ko,v are constants, substitution from (|4.8|) into H4.9|l yields the condition 



(4.10) 



2aa - a-" - 4w^ + 4 ^o + — = 



This nonlinear equation can be exactly solved in terms of the series expansion 
(see the Appendix) 



2k+l 



(4.11) 



fc=0 



where all the coefficients of even index, a2k , are zero and 

[2ik + l)-j] 



2fe+l 

a2,+i=(-i)v n 



{2k + 3-j)[j{j-2) + X^] 



ai 



where 



0.1 



2wo 



Provided that > 1, we thus get 



m^mo'^ a2fe+i t 

k=0 



2fc+l 



= — j} a2kt 
ait ^-^ 

k=0 



2k 



(4.12) 



(4.13) 



(4.14) 



where (see e.g. [TBj 'l 



-1 
ai 






ai 














03 





ai 














as 





ai 








as 





as 









Ok+l 


flfe 












It is worthwhile noting that for >> jt^ one would get from H4.10|l the 



solution 



a = ci + W cf 



"o "-0 



sm 



41^2 fc2 (i + c2) 



where ci and C2 are constants. Inserting such a solution into (|4.8|l . H4.1U|I . and 
solving the equation of motion (|1.2II we would obtain 



q(t) = C\ sin • 



4tan"^ 



"202 IP- 



1 - Cli 



tan \l^lkl (i + C2) 
(4.15) 



where Ci and (/)o are real constants. Due to H3.29|l . the phase Q now reads 

t 



e = 2wo 



fc=l 



(4.16) 



In the special case in which /i = (fco = 0), we simply have 

a{t) ^ait . (4.17) 

Inserting (|4.17|) in (|4.1UI) provides just the condition (|4.13|l . while the equation 
of motion p.2|l becomes 



2 ,,2 



i\ 1 



4/ i2 



9 = 



(4.18) 
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which affords the generai soiution 



q = C,t''+ +C2t^- , /3± =±yi-f^2^2 . (419) 

At this stage some comments are in order. In doing so, it turns out to be 
meaninfui to go back; to the originai form of the auxihary equation H1.6|l . which 
now reads 

1 



a + nl i^k^, + ^ja~ul^ = . (4.20) 

Next, we consider the equation of motion corresponding to H4.8|l . That is (see 

q+-q + uj\t)q^ , (4.21) 
a 

where w^(t) is given by 

due to p.5|l and (|4.9|) . By means of the change of variable q — y / y/a Eq. H4.21|l 
is transformed into the equation 

2/ + f^offco + ^) y - , (4.23) 
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which can be led to a Bessel function. For example, if ^ = ^oko, H.qL'^ = \ ^ P' 
Equation H4.23|l is solved by 

y = VtZp{tt) , (4.24) 

where Zp{t) denotes a Bessel function of order p defined by the differential 
equation |16j 



We notice that when the quantity is negligible with respect to the other 
terms, Eq. (|On|) can be approximated by an equation of the type (14.23(1 . 
namely 



a + % \ H + -pT] ^^ = . (4.26) 



5 Discussion 



We have formulated the uncertainty product (A^Q) (A^P) for the generalized 
(time-dependent) oscillators (TDO's) with variable mass and frequency, where 
Aq,(5 and AqP are the variances between extended coherent states \a) of the 
conjugate position and momentum operators Q and P. The basic starting point 
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of our study is constituted by the auxiliary equation l|1.6f) which is associated 
with the equation of motion of the TDO under consideration by formula 1)1.71) . In 
the case of the conventional harmonic oscillator, we have simply found that Eq. 
1)2. 7|l gives rise to two (real) exact different solutions, (|2.8|) and (|2.9|) . whose the 
first minimizes exactly the uncertainty expression H3.12|l (M — 0, & — 0), while 
the second does not minimize exactly this expression, but only approximately. 

Another case in which the auxiliary equation H1.6|l can be explicitly solved 
is represented by the particular Kanai-Caldirola oscillator l|2.12|l with 57^ (t) = 
f2§ = Wq ^ 7^/4 = const and m(t) = TOoexp(7t) (7 ~ Mq is the damping). 
We obtain two possible exact solutions of the corresponding auxiliary equation 
accordingly to fig > or fig < 0. One solution (for fig > 0) is given by 
a = , ^ and, therefore, it resembles what happens for the usual harmonic 

oscillator (see (|2.8ll '). while the other (for Jig < 0) is furnished by (|2.15|) and 
does yield minimum states only approximately. 

As one expects, the uncertainty relation H3.22|l can even be written as 



V(0|Q2|o) (0|P2|0) = - 



9 1 M ^ ^ 
1 + 4cr (7 - — cr 



(5.1) 



2 

where the expectation values between the vacuum state of and are 

(0|Q2|o) = la2 , {0\P^\0) = hm\^\' , (5.2) 
m 

with 1^1^ expressed by ()3.9|) . Since in the case of minimum uncertainty a = 
cm}/'^ (see (|3.24ll '). the expectation values (0| |0) and (0| |0) turn out to 
be constant, namely 

h 

A? 

so that 



(0|g2|0) = ftc2 , (0|p2|o)= , (5.3) 



h 



v/(0|g2|0) (0|P2|0) = - (5.4) 

corresponds to the minimum uncertainty. In this situation, from H3.16|l we infer 
that the Bogolubov coefficients /i(t) and v{t) are identical to those characterizing 
the usual harmonic oscillator: 

^i{t) ^ 1 , v{t) ^ . (5.5) 

On the other hand, in the light of the criterium (|3.27|l the Hamiltonian operator 
H3.1|l reads 

H{t)^^H{0) (5.6) 
m(t) 

where muj = moojo ~ l/(2c^). We remark that the vacuum state carries of 
course zero momentum (0| P |0) =0 but not zero energy, in the sense that the 
requirement (|3.27|) entails 

(0| Hit) |0) ^ (0| H{0) 10) = ^ , (5.7) 
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which can be interpreted as the zero-point energy of the time-dependent quan- 
tum oscihator minimizing uncertainty states. 



Appendix 

In order to solve the nonUnear equation 



2aa-a^ -410^+ i^H^ + —j =0 (A.l) 
the series expansion 

oo 

(A.2) 

k=0 

can be introduced. From l|A.2|l we thus obtain 



J^ht" , a'' = J2^kt'' , aa = J2dkt'' (A.3) 

k=l k=0 k=0 

where 

k k 

&fc = ^(j + 1) (fc - j + l)aj+i afc_j+i , Ck = '^ajak-j (A. 4) 

J=0 3=0 

and 

k 

dk=^[2 + k- j){l + k- j)ajak+2-o ■ (^-5) 
i=o 

Inserting (|A.3|I into ljA.l|l and coUecting terms of same power of t yields 

A^co = A^ci = , 

2do - 6o - 4wg + M^cq + C2 = , 

2d„ - hn + C„ + A^ Cn+2 = 

Vn > 1. In terms of the a^'s we have 
A^ ao^ = 2 A^ ao ai = , 

4 flo 02 — fli^ — 4 cjo^ + flo^ /^^ + A^ (ai^ + 2 ao = , 
12 Oo 03 + 2 /i^ ao ai + 2 A^ (ai 02+ ao as) = , 

6 ai as + 24 ao + j.P' (ai^ + 2 ao a2) + A'^( a2^ -f 2 ai as -I- 2 ao 04) = , 
16 ai a4 + 4 02 as + 40 ao as -|- 2 (ao as + ai a2)-l- 
+ 2A^ (ai a4 -|- 02 as + ao as) = , 
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30 oi 05 + 12 02 04 + 60 oo 06 + 3 03^ + /i^ (2ao 04 + 2 oi 03 + 02^)+ 

+ (03^ + 2 01 0,5 + +2 02 04 + 2 00 fle) = , 
48 Oi og + 24 02 05 + 12 0,3 0,4 + 84 Oo 07 + 2 fj,^ (oq 05 + Oi 04 + 02 03) + 

+ 2 oi og + 02 05 + 03 04 + Oo 07) = , 
70 oi 07 + 40 02 oe + 22 03 05 + 112 oq os + 8 04^ + /x^ (03^+2 oo og + 2 oi 05+ 

+ 2 02 04) + (04^ + 2 oi 07 + 2 02 06+2 03 05 + 2 oq os) = , 
24 04 05 + 96 oi og + 60 02 07 + 36 03 oe + 144 oq 09+^^(2 oq 07 + 2 oi 06+ 

+ 2 02 05 + 2 03 04) + A^(2 04 05 + 2 oi og + 2 02 07+2 03 oe + 2ao 09) = , 
36 04 06 + 126 oi og + 84 02 og + 54 03 07 + 180 oq oio + 15 05^+ 

+ /(X^ (04^+2 Oq Og + 2 Oi O7 + +2 02 06 + 2 O3 05) + 

+ A^( 05^ + 2 04 06 + 2 Oi 09 + 2 02 Og + 2 03 07 + 2 oo oio) = , 
52 04 07 + 40 05 Og + 160 Oi Oio + 112 02 09 + 76 03 og + 220 Oq aii + 
+ fi^{2 Oq 09 + 2 04 05 + 2 oi Og + 2 02 07 + 2 03 05)+ 

+ A^(2 O4 O7 + 2 O5 Og + 2 Oi Oio + 2 02 Og + 2 03 Og + 2 Oq aii)+ — , 

72 04 Og + 54 05 07 + 198 Oi On + 144 02 oio + 102 03 09 + 264 oq 012 + 24 06^+ 
+ u'^ 06^ + 2 /i^ (05^ + 04 O6 + Oo Oio + oi 09 + 02 Og + 03 07) + 
+ 2A^( 04 Og + 05 07 + oi Oil + 0,2 aw + O3 Og + Oo 012) = 

etc. Since cq = Og then we must require oq = 0, which also imphes Ci = 
2ooOi = 0. Moreover, from the third equation we obtain oi (7^ 0), that is 
Oi^(A^ — 1) = 4cjo^- Setting oo = drastically simplify the whole system in 
that, due to the typical bilinear structures involved, the net consequence is that 
coefficients Ofe of the even index must vanish as well. For instance, once we 
set Oq = the fourth equation above reads 2 A^ 0102 = 0. Further insertion of 
02 = into the system results into the sixth equation (A^ + 8)0104 = and so on. 
As regards the series coefficients of the odd index, they can be determined 
recursively from the the fifth, the seventh etc., equations of the system. These 
equations now read 

oi {fi^ oi + 6 03 + 2 j/^ 03) = , 

30 Oi 05 + 3 03^ + 2 /K^ oi 03 + A^ (03^ + 2 ai 05) = , 

70 oi 07 + 22 03 05 + /U^ (03^ + 2 oi 05) + 2 A^( oi 07 + 03 05) = , 

126 oi 09 + 5403 07 + 15 05^ + 2/i^ (oi 07 + 03 05) + A^(a5^ + 2 oi ag+ 

+ 2 03 07) = , 
5405 07 + 198 oi Oil + 1020309 + /x^(2 oi 09 + 2 03 07 + 05^)+ 

+ 2A^( 05 07 + oi On + 03 09) = , 
94 05 Og + 286 oi 013 + 166 03 on + 35 07^ + 2 /n^ (oi On + 03 09 + 05 07)+ 

+ A^ (07^ + 2 a5 09 + 2 oi ai3 + 2 a3 On) = , 
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etc. It is an easy matter to handle this system iteratively, thus giving all the 
odd index coefficients in terms of the ai. Precisely, a httle bit of inspection 
reveals the existence of the recursion relation 14.12|l . which can be subsumed 
even into 
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